Measures of quantum synchronization in continuous variable systems 
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We introduce and characterize two different measures which quantify the level of synchronization 
of interacting continuous variable quantum systems. The two measures allow to extend to the 
quantum domain the notions of complete and phase synchronization. The Heisenberg principle sets 
a universal bound to complete synchronization. The measure of phase synchronization is in principle 
unbounded, however in the absence of quantum resources {e.g. squeezing) the synchronization level is 
bounded below a certain threshold. We elucidate some interesting connections between entanglement 
and synchronization and, finally, discuss an application based on quantum opto-mechanical systems. 



In the 17th eentury, C. Huygens noticed that the oscil- 
lations of two pendulum clocks with a common support 
tend to synchronize (Fig. flla) jlj. Since then, analo- 
gous phenomena have been observed in a large variety 
of different contexts, e.g. neuron networks, chemical re- 
actions, heart cells, fireflies, etc. |H|. They are all in- 
stances of what it is called the spontaneous synchroniza- 
tion effect where two or more systems, in the complete ab- 
sence of any external time-dependent driving force, tend 
to synchronize their motion solely due to their mutual 
coupling. The emergence of spontaneous synchroniza- 
tion in so many different physical settings encouraged its 
investigation within classical non-linear dynamical sys- 
tems. Here, given the time evolution of two dynami- 
cal variables, like the position of two pendula, standard 
methods exist to verify whether their motion is synchro- 
nized [2]. For quantum systems however the same ap- 
proaches cannot be straightforwardly extended due to 
the absence of a clear notion of phase space trajectories. 
The aim of this work is to address this problem, develop- 
ing a consistent and quantitative theory of synchroniza- 
tion for continuous variable (CV) systems evolving in the 
quantum regime [3] . To this aim we introduce two differ- 
ent quantum measures of synchronization extrapolating 
them from notions of complete and phase synchroniza- 
tion introduced for classical models. We will show that 
quantum mechanics set bounds on the achievable level 
of synchronization between two CV systems and we will 
discuss the relationship between entanglement and syn- 
chronization. We finally apply our approach for studying 
the dynamics of coupled opto-mechanical systems [H [S] . 

In the quantum domain synchronization has been stud- 
ied in various contexts, like quantum information proto- 
cols [6], two- level systems [7] and stochastic systems [8]. 
While our measures could in principle be extended also to 
these cases, our endeavor is specifically framed in the re- 
search line investigating the spontaneous synchronization 
of micro- and nano-mechanical systems [5HT5]. Recent 
experimental advances allow to realize opto-mechanical 
arrays composed of two or more coupled mechanical res- 
onators controlled close to their quantum regime by laser 
driving [TBHin]- Such devices have all the properties 




FIG. 1. (Color online) Original Huygens' sketch [Tl of two syn- 
chronizing pendulum clocks (a) and the quantum mechanical 
analogue consisting of two (or more) coupled opto-mechanical 
systems (b). Here, mechanical resonators are driven into 
self-sustained oscillations by the non-linear radiation pressure 
force of independent optical modes. A weak mechanical in- 
teraction is responsible for the spontaneous synchronization 
of the limit cycles. All symbols are defined in the main text. 



(non-linear dynamics, limit cycles, etc.) which are nec- 
essary for the emergence of spontaneous synchronization 
[H [20] and indeed some first experimental evidences of 
this effect have been found [HI [15] . 



Quantum synchronization measures:- In a purely clas- 
sical setting, synchronization is mostly studied in the 
context of autonomous non-linear systems undergoing 
limit cycles or chaotic evolution (linear systems being 
usually excluded because they converge to constant or 
unstable solutions). In this scenario one can identify 
different forms of synchronization j2]. Complete syn- 
chronization is achieved when (say) two subsystems Si 
and ^2, initialized into independent configurations, ac- 
quire identical trajectories under the effects of mutual 
interactions. Specifically, given two CV classical sys- 
tems characterized by the (dimensionless) canonical vari- 
ables qi (t) , pi (i) and 52 (t) , P2 (t) describing the evolution 
of Si and 6*2 in phase space, complete synchronization is 
reached when the quantities q-{t) :— [qi{t) — q2{t)]/^ 
and p-{t) :— [pi(i)—p2(^)]/'\/2 asymptotically vanish for 
large enough times |21j . Phase synchronization is instead 
achieved when, under the same conditions detailed above, 
only the phases <fj{t) = a.Yctan[p j{t)/qj{t)] are locked: 



i.e. when the quantity ip-{t) := ^i{t) — 'P2{t) asymptot- 
ically converges to a constant phase shift ipo & [0, 2n]. 

One can already foresee that extending the above con- 
cepts to quantum mechanical systems is not straight- 
forward and that some fundamental limits could exist 
that prevent the exact fulfillment of the conditions given 
above. In particular, identifying the dimensionless quan- 
tities qj{t), Pj(t) as quadrature operators obeying the 



canonical commutation rules [qj{t),pj'{t)] 
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the relative coordinates q- (t) and p_ (t) will correspond 
to generalized position and momentum operators of the 
same (anti-symmetric) mode of the system. Accordingly 
the uncertainty principle will now prevent the possibility 
of exactly achieving the condition required by classical 
complete synchronization. 

To turn this into a quantitative statement, we identify 
g_ (i) and p- (t) as synchronization errors and introduce 
the following figure of merit 



ss)--={q-{tr+p-itr) 
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(1) 



gauging the level of quantum complete synchronization 
attained by the system (here (■ ■ • ) implies taking the ex- 
pectation value with respect to the density matrix of the 
quantum system). We then observe that the Heisenberg 
principle requires ((7_(f)^)(p_(t)^) > 1/4 and hence 



Sc{t) < 



1 
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which sets a universal limit to the complete synchro- 
nization two CV systems can reach. On the contrary, 
in a purely classical theory, Sc{t) is in principle un- 
bounded [55]. A small value of Sc{t) can have two pos- 
sible origins: the mean values of q- (t) and p- (t) are not 
exactly zero, and/or the variances of such operators are 
large. The former situation can be interpreted as a clas- 
sical systematic error (the trajectories of the limit cycles 
are synchronized but not exactly equal) while the latter is 
due to the influence of thermal and quantum noise. The 
classical systematic error can be easily excluded from the 
measure of synchronization by using the same expression 
of Eq. (fTl) but after the application of the change of vari- 
ables 

Q_ (i) ^ g_ (t) - (g_ (t)) , p_ (t) ^ p_ it) - (p_ (t)} . (3) 

This gives a relative measure of synchronization which is 
always larger than the previous absolute one and which 
may be preferable whenever the aim is that of selectively 
investigating purely quantum mechanical effects. Obvi- 
ously, the bound of Eq. pi holds also for this relative 
measure. 

Constructing a quantum analogue of the phase syn- 
chronization condition is more demanding due to the con- 
troversial nature of the quantum phase operator(s), see 
e.g. Ref. [53]. We adopt a pragmatic point of view that 
allows us to target departures from the ideal (classical) 



condition, due to quantum fluctuations. To do so, we 
write the operator aj{t) := [qj{t) + ipj{t)]/\/2 of the j-th 
system as 



a,(t) = h(t)4-4(i)]e^'^^-W, 



(4) 



where rj{t) and (pj{t) are the amplitude and phase of 
the expectation value of aj{t), i.e. {aj{t)) = rj(i)e*'^j(*^ 
With this choice, the Hermitian and anti-Hermitian part 
of a'j{t) — [qj{t) + ip'j{t)]/^/2 can now be interpreted 
as fluctuations of the amplitude and of the phase re- 
spectively (indeed this is the reason why in quantum 
optics q'j{t) and p'j{t) are often called amplitude and 
phase quadratures). If two CV systems are on aver- 
age synchronized such that the phases of {ai(t)) and of 
{a2{t)) are locked, then the phase shift with respect to 
this locking condition can be associated to the operator 
p'_{t) — [p'i{t) — P2(^)]/v2. A measure of quantum phase 
synchronization can then be obtained through the quan- 
tity 



S,{t) := -{p'_{tn 



(5) 



Differently from the measure (fTl), Sp can be in principle 
arbitrarily large. Nonetheless, even in this case, quantum 
mechanics sets a non trivial bound relating the quantum 
optical notion of non-classicality with phase synchroniza- 
tion. If two CV quantum systems evolve in time such 
that their P-function [S] [24] is always positive, then per- 
fect phase synchronization is impossible and one has 



positive P- function => Sp{t) < 1 



(6) 



Indeed a value of (p'_(i)^) below 1/2 implies the existence 
of collective squeezing and so the impossibility of a phase 
space representation of the state through a positive P- 
function. 

Furthermore the specific structure of the limit cycles 
associated with the average quantities rj{t) and ipj{t) 
may lead to additional bounds for Sp. If, for example, 
i) the system under consideration exhibits mean values 
quantities {aj{t)) converging to approximately circular 
limit cycles in the phase space, ii) the noise operating 
in the system is not phase sensitive {i.e. is invariant for 
phase space rotations) and Hi) the interaction potential 
between the two systems is of the form Hint — — /i(aia2 + 
a2a\); then it is reasonable to conjecture that (p'_{t)'^) > 
{q'-{t)'^). This together with the Heisenberg principle, 
leads to the bound 



Sp{t) < S,{t) < 1 



(7) 



While referring to the Supplementary Material for an 
heuristic derivation of Eq. ([7|, we remark that such in- 
equality is consistent with the results shown later on 
opto-mechanical systems. 

Quantum correlations and synchronization: - Synchro- 
nization and entanglement are both associated with the 



presence of correlations between two or more systems. It 
is thus natural to ask if, in the quantum regime, there 
is a strong interplay between the two effects. Quite sur- 
prisingly however it turns out that, according to our mea- 
sures, the stationary state of two CV systems can possess 
maximum amount of complete or phase synchronization 
without being necessarily entangled. For instance a sys- 
tem converging to two factorized coherent states evolving 
in time such that (ai(t)) = {a2{t)), exhibits maximum 
complete synchronization (Sc — 1) but has no entan- 
glement. Similarly consider two locally squeezed states 
rotating in phase space such that (ai(i)) — (02 (i)) and 
(pi(O^) = {Pii^Y) = €j with p'j, being the quadrature 
orthogonal to the phase space cycle of subsystem k as 
defined in Eq. Q (said in simpler words, these are two 
squeezed states moving like synchronized clock hands 
in phase space). This state has arbitrary high phase 
synchronization Sp = \^~^ but it is clearly not entan- 
gled. Entanglement appears hence to enforce correlations 
which are qualitatively different from those required to 
yield high values for Sc{t) and Sp{t). A better insight 
on this can be obtained by considering the very precur- 
sor of all CV entangled states, i.e. the EPR state [21] 
which describes the ideal scenario of two systems having 
same positions but opposite momenta. It is thus clear 
that synchronization requires different constraints which 
could have instead a relationship with other measures of 
quantum correlations like quantum discord (see e.g. our 
successive results on opto-mechanical systems). We con- 
clude this section with an open question on the converse 
problem: what kind of synchronization phenomenon cor- 
responds in the quantum limit to EPR correlations? EPR 
entanglement could be identified as a mixture of com- 
plete and anti-synchronization, i.e. qi(t) = 92 (i) and 
Pi{t) = —p2{t). Recently this unconventional regime 
called mixed synchronization has been introduced and ob- 
served in classical non-linear systems 30], but whether 
this concept is relevant and extendible in the quantum 
domain is still unexplored. 

Measures and bounds at work:- Opto-mechanical de- 
vices [U [S] provide the perfect setting where our mea- 
sures for synchronization can be directly applied. We 
thus identify 5"! and S2 with two approximately identi- 
cal mechanical resonators (see Fig. [lib) coupled to in- 
dependent cavity optical modes (needed to induce self- 
sustained limit cycles) and mutually interacting through 
a phonon tunneling term ^ of intensity p,: 



H = Y^ [Aja]aj + Ujh]hj - ga]aj{bj + b]) 



iE{aj - a])] - m(6i&^ + fe^6i) 



Kt; 



(8) 



In this expression, for j — 1,2, aj and bj are the op- 
tical and mechanical annihilation operators, tUj are the 
mechanical frequencies, Aj are the optical detunings, g 
is the opto-mechanical coupling constant, while E is the 



laser intensity which drives the optical cavities {h = 1). 
For simplicity g and E are assumed to be equal in both 
systems while wi and uj2 can be slightly different. Dissi- 
pative effects are included adopting the Heisenberg pic- 
ture and writing the following quantum Langevin equa- 
tions ^E\, 

hj = [-K - iAj + ig{bj + b])]aj + E + -j2Ka"\ (9) 
bj — [—7 - iuij\bj + iga}^aj -I- ifib^^j + y^b™'. 

Here k and 7 are, respectively, the optical and me- 
chanical damping rates while a*" and 6"' arc the in- 
put bath operators. These are assumed to be white 
Gaussian fields obeying standard correlation relations, 

(a;"(i)ta™(i') + a™(i')a"(0^) = %"5(t - t') and 
(6j"(i)tK.?(t') + 6;-?(i')6!"(i)t) = (2nfc + l)S,y6{t ~ t'), 



where rit, — [exp( 






1] ^ is the mean occupation 



number of the mechanical baths which gauges the tem- 
perature T of the system [5S] (since we are only interested 
in the situation in which uji ~ u;2, the parameter rif, can 
be safely taken to be equal for both oscillators). 

In the semi-classical regime |25| we express each an- 
nihilation or creation operator 0{t) of (|9| as sums of 
mean values contributions {0(t)) (representing the clas- 
sical trajectories of the system) plus small fluctuations 
terms 0'{t), i.e. we write 0{t) = {0{t)) + 0'{t). This 
allows us to determine, in a semi-classical approxima- 
tion, the expectation values {0(t)) in terms of a set of 
classical non-linear differential equations [5S]. Setting 
Aj = —ujj (driving detuning) and choosing the laser am- 
plitude E of Eq. Q large enough, we make sure that such 
solutions yield limit cycles as classical steady state con- 
figurations (for more details on opto-mechanical multi- 
stability and limit cycles see |26j ) . Quantum fluctuations 
are instead obtained by computing the covariance ma- 
trix C(i), with entries given by Ci^t{t) = {Ri{t)Ri{t)^ + 
Re{ty Ri{t))/2, the expectation value being taken on the 
initial state and Ri being the components of the vector 
R = {a[, a'l , b'l , b'^ , a'2 , a'2 , &2 1 ^2 ) ■ I^ particular this gives 
us direct access to the mechanical variances {q-{t)'^) and 
(p_(t)^) which define the complete synchronization level 
via Eq.Q. By applying the linearization procedure, we 
implicitly performed the change of variables correspond- 
ing to Eq. ([3]) and so we automatically excluded the sys- 
tematic synchronization error due to slightly different 
average trajectories. As a consequence the only source 
of disturbance bounding our measure of synchronization 
will be quantum (or thermal) fluctuations. 

Estimating phase synchronization as in Eq. ([5]) requires 
instead a further step as the latter has been defined with 
respect to a reference frame rotating with the phases of 
the average trajectories, see Eq. Q. This corresponds 
to a diagonal and unitary operation on R, built up on 
the phases ipaAt) = arg(oi(i)), (/^ajlO = arg(a2(i)), 
etc., of the classical orbits: i.e. R —>■ R' = U{t)R with 
U{t) = diag[e-''^"i(*\e*'^°i (*),•••]. The associated co- 
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FIG. 2. (Color online) Subfigure (a): simulation of the com- 
plete (blue) and phase (green) synchronization measures M 
and ([5| between the mechanical resonators as functions of 
time (m units of r = 2-k/ijJi). The dashed lines indicate the 
corresponding time averaged asymptotic values, i.e. the quan- 
tities Sx ~ limT-»oo ^ Jg Sxit)dt for x = c,p. Setting tJi = 1 
as a reference unit of frequency, the other physical parame- 
ters which have been used in the simulation are: cj2 = 1.005, 
7 = 0.005, Aj = -cjj, K = 0.15, g = 0.005, /i = 0.02, m = 
and E = 320. Subfigure (b): time averaged complete (cir- 
cles) and phase (squares) synchronization and Gaussian dis- 
cord Dg (diamonds) as functions of the coupling constant 
fj,. Subfigure (c): time averaged synchronization measures as 
functions of the bath mean phonon number nt. Subfigure (d): 
Synchronization between two arbitrary mechanical modes of 
a chain of 20 coupled opto-mechanical systems as a function 
of the lattice distance h. All subsystems are assumed to have 
the same mechanical frequency lj = 1. 



variance matrix is C'{t) = U {t)C {t)U {t)\ from which 
we can directly extract the mechanical variance {p'^{t)) 
entering Eq.([5]). 

A simulation of the complete and phase synchroniza- 
tion between the mechanical modes is plotted in Fig. 
[2Ja using realistic values for the parameters [H |S] (see 
caption for details). After an initial transient, the sys- 
tem reaches a periodic steady state in which Sc{t) and 
Sp{t) are significantly larger then zero, implying that 
both complete and phase synchronization take place in 
the system. Their value is consistent with the funda- 
mental limit (pi) imposed by the Heisenberg principle 
and with the heuristic bound ([7]) presented in the pre- 
vious section. Indeed we numerically find that quantum 
squeezing in the p'_{t) quadrature, needed to overcome 
the non-classicality threshold ^ , is absent in the system. 
Fig. [2]b and Fig. [2jc report instead the behavior of the 
time averaged measures of complete and phase synchro- 
nization for different values of the coupling constant and 
of the bath temperature. We vary n from zero [57] to a 
maximum threshold above which the classical equations 
are perturbed too much destroying the limit cycles. 



Finally we have checked if quantum correlations are 
present in the system verifying that, consistently with 
the difference between entanglement and synchronization 
detailed in the previous section, for many choices of the 
parameters entanglement negativity is always zero even 
though synchronization is reached. On the contrary, non- 
zero level of Gaussian quantum discord [23] (Fig. [2]b) 
between the two mechanical modes is observed for all 
values of /i that lead to synchronization. Still our data 
are not sufScient to clarify the functional relationship 
between discord and synchronization (if exists) . 

The synchronization observed between the oscillators 
is expected to emerge also when more than two parties 
are present in the setup. In particular we focus on the 
case of a (closed) chain formed by N opto-mechanical 
systems with first neighbor interactions (the Hamiltonian 
being the natural generalization of ([8]) with uniform pa- 
rameters). As before, we enforce the driving detuning 
condition A = -co and set the laser intensities E in or- 
der that each opto-mechanical system converges to a sta- 
ble limit-cycle. Once these prerequisites are fulfilled we 
linearize the dynamics around the classical steady state, 
which is assumed to be the same (synchronized) in each 
site. I.e. {aj{t)) = a(t) and {bj(t)) = I3{t) for aU j. This 
corresponds to a mean-field approximation applied only 
to the classical dynamics, while the fluctuation terms a' 
and 6'- can be treated exactly (without mean-field) since 
the associated Hamiltonian is quadratic. Fig.[2jd reports 
the results obtained for two mechanical modes separated 
by h lattice steps: we notice that the synchronization 
level among the various elements persists even if an ex- 
ponential decay in h is present (a behavior which is con- 
sistent with the one-dimensional topology induced by the 
selected interactions) . 

Summary:- We have quantitatively studied the phe- 
nomenon of spontaneous synchronization in the setting 
of coupled CV quantum systems. We have shown that 
quantum mechanics sets universal limits to the level of 
synchronization and discussed the relationship between 
this phenomenon and the emergence of quantum corre- 
lations. Finally we have analyzed the spontaneous syn- 
chronization of opto-mechanical arrays driven into self- 
sustained oscillations. A large number of open aspects 
are worth being further investigated, among which: the 
interplay between quantum correlations and synchroniza- 
tion, the application of this theory to other physical sys- 
tems and the interpretation of synchronization as a useful 
resource for quantum communication and quantum con- 
trol. 
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SUPPLEMENTAL MATERIAL 

The hypothesis underlying the conjecture Q are: (i) 
the system admits hmit cycles which are approximately 
circular {aj{t)) ~ rj-e"*"^*, (ii) thermal (or quantum) 
noise is not phase sensitive in the sense that it is in- 
variant for phase space rotations, and (iii) the interac- 
tion potential between the two systems is of the form 
Hint = — /i(aia2 -l- a2ai). These assumptions are often 
valid for optical or mehcanical modes under the rotating 
wave approximation. For such systems the interaction 
can be written in terms of symmetric and anti-symmetric 
normal modes Hint ~ —fia_^a-\- + fioLa- which up to 
a renormalization of the bare frequencies (oj — > cj — /i) 
is equivalent to -ffj'„t = 2/ia_a_. This means that, be- 
cause of the interaction, the anti-symmetric mode rotates 
with a frequency 2/i faster with respect to the symmetric 
mode. 

Since the phase of each individual limit cycle is arbi- 
trary, the same must be true for the symmetric mode. 
Unless anti-synchronization appears, one can guess the 
following structure for the linearized equations for the 
fluctuations: 



di 






-leff 





1+ 



(10) 



The reason is that the symmetric mode must be stable 
in the amplitude (with effective damping 7e// > 0) but 
its phase should be completely free to diffuse (the corre- 
sponding Lyapunov exponent must be zero). 

The dynamics of the anti-symmetric mode, linearized 



around a synchronized solution, will be like Eq. ( 10 ) plus 



a frequency shift of 2/j, due to the interaction potential: 



d 
dt 



-leff 

2n 



-2fi 




(11) 



The new matrix, for /i 7^ 0, has negative eigenvalues and 
this fact is the origin of synchronization. Since we as- 
sumed the noise to be phase insensitive the diffusion ma- 
trix must be proportional to the identity. In this case 
one can easily find the steady state and check that in- 
deed (p_?) > (q^), and hence the bound i7h holds. 

Of course this heuristic argument is very hand-waving 
but it gives the physical intuition that, in a classical or 
quantum system, the precision of phase synchronization 
may be bounded by the precision of amplitude synchro- 
nization. 



